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Investigation  of  Electromagnetic  Field  Propagation  in 
Reverse  Saturable  Absorbers 

1.  Introduction 

The  general  area  of  nonlinear  absorbers/optical  limiting  has  grown  substantially 
over  the  last  two  decades.  There  is  increasing  interest  in  various  mechanisms  of  optical 
limiting.  A  number  of  review  articles  and  conference  reports  demonstrate  the  vitality  of 
this  subject  [1]-[11].  In  general,  an  optical  limiter  keeps  the  power,  intensity,  energy  or 
energy  density  transmitted  by  an  optical  system  below  a  predetermined  maximum  value 
that  is  independent  of  the  size  of  the  input  pulse  while  maintaining  a  high  transmittance  at 
low  input  power.  The  many  applications  of  the  device  include  laser  power  regulation, 
laser  mode-locking,  optical  pulse  shaping,  signal  level  processing,  and  sensor/detector 
protection. 

Importantly,  the  device  must  operate  within  a  fast  enough  response  time  to  protect  a 
sensor  or  regulate  pulsed  laser  sources.  The  transmission  must  be  reduced  in  a  time 
period  shorter  than  the  pulse  width.  With  laser  pulse  durations  ranging  from  milliseconds 
to  femtoseconds,  electrooptic  devices  such  as  electrooptic  switches  that  operate  on  the 
nanosecond  time  scale  are  often  too  slow.  On  the  other  hand,  passive  devices  use  a 
nonlinear  optical  medium  that  has  a  fast  response  time  and  acts  as  a  sensor,  processor  and 
modulator  [8],  However,  typical  nonlinear  materials  are  not  sufficient  optical  limiters  by 
themselves.  As  a  result,  the  passive  system  places  the  nonlinear  optical  components  in  or 
near  a  focal  plane.  The  focal  point  should  be  located  near  the  exit  plane  of  the  material  to 
utilize  the  entire  path  length. 

The  laser  threat  occurs  at  wavelengths  within  the  transparency  window  of  the 
atmosphere.  Thus  the  most  dangerous  laser  sources  will  be  in  the  following  windows: 
visible-near  infrared  region  between  0.4|im  and  1.4pm;  infrared  region  I  (IR  I)  window 
between  1.5pm  and  2.7pm,  infrared  region  II  (IR  II)  window  between  3pm  and  5.8pm 
and  infrared  region  in  (IR  III)  window  between  7.5pm  and  14pm.  The  optical  limiters 
working  inside  these  wavelength  ranges  are  of  unparalleled  important.  The  main  types  of 
detectors  fall  into  two  categories;  namely,  direct  (eye)  and  indirect  (photo  detector,  image 
intensifiers  and  IR  night  vision  sensors). 

The  maximum  permissible  transmitted  energy  is  dependent  on  many  factors 
including  the  optical  wavelength,  pulse  width,  linear  transmittance,  maximum  output, 
frnumber  of  the  system,  damage  threshold  and  stability.  The  various  situations  not  only 
cover  a  very  wide  range  of  wavelengths,  which  includes  infrared  and  the  whole  visual 
range  but  also  require  different  types  of  optical  limiter  mechanisms,  e.g.  fluence  limiter, 
energy  limiter  and  irradiance  limiter. 

Optical  limiting  now  requires  a  sophisticated  set  of  properties  including  various 
nonlinear  materials,  different  physical  processes,  various  device  designs  including  a 
number  of  lenses  and  other  optical  devices,  optical  pulse  shaping,  and  the  characteristics 
of  the  laser  source  including  wavelength,  pulse  duration  and  pulse  repetition  rate.  This 
range  of  materials,  processes,  and  devices  involve  complex  optical  propagation  requiring 
a  vast  diversity  of  parameters  that  cannot  be  fully  analyzed  without  the  use  of  various 
numerical  codes.  Along  with  new  materials  and  device  structures,  new  numerical  codes 
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will  be  key  elements  in  the  development  of  optical  limiters.  Numerical  investigation  will 
be  required  to  make  additional  advances  by  examining  the  various  combinations  of 
material,  physical,  and  device  parameters  that  would  be  too  time  consuming  and/or  costly 
to  perform  in  the  laboratory. 

Physical  mechanisms  used  in  optical  limiters  include  nonlinear  refraction,  nonlinear 
scattering,  two  photon  absorption,  thermal  nonlinearity,  and  reverse  saturable  absorption. 
Materials  used  in  optical  limiters  include  fullerenes,  porphyrins,  liquid  crystals,  carbon 
black  and  nanostructured  materials  [1]-[10],  Various  device  design  and  structures  also 
play  an  important  role.  Examples  of  different  device  design  include  tandem  and  graded 
density  structures  [7], 

One  of  the  first  passive  limiter  was  developed  by  Leite  et.  al.  [12]  in  1967  based  on 
the  thermal  lensing  effect  in  nitrobenzene.  This  device  regulated  the  output  power  of  a 
continuous  wave  (CW)  Argon  ion  laser.  It  wasn’t  until  1983  that  CS2  [13]  was  used  as  an 
optical  limiting  material  for  nanosecond  pulses.  The  experiment  used  a  focused 
geometry  and  showed  that  the  mechanisms  involved  in  limiting  transmission  were 
catastrophic  self-focusing  together  with  absorption  and  scattering  resulting  from  laser 
induced  breakdown.  In  the  last  two  decades  experiments  and  theory  on  optical  limiting 
has  greatly  expanded  with  the  mechanisms  including  two-photon  absorption,  excited  state 
absorption,  nonlinear  refraction  and  nonlinear  scattering.  A  variety  of  materials  are  being 
investigated  and  characterized  to  improve  the  performance  of  optical  limiters.  Recent 
research  interests  are  focused  on  organic  material,  like  dyes,  liquid  crystals  and 
nanostructures,  which  exhibit  very  high  nonlinearity. 

The  damage  threshold  is  determined  by  several  factors  including  the  total 
energy/fluence  entering  the  eye  [14]-[19],  the  pulse  duration  [14],  the  spot  size  [16],  and 
the  number  of  incident  pulses  [15].  Ocular  media  are  transparent  in  the  range  between 
0.35pm  to  1.3pm  [14]  Early  experiments  by  Birngruber  et.  al.[17]  measured  laser  eye 
damage  at  the  ED50  level  (the  50%  probability  of  retinal  lesion  development)  as  a 
function  of  pulse  width  from  about  0. 1  sec  to  about  100  fs.  For  pulse  durations  above  the 
microsecond  level  (~20  ps),  the  damage  mechanism  is  due  to  thermal  effects  [14]  and  is 
irradiance  dependent  (i.e.  W/cm2).  In  the  pulse  duration  range  from  nanoseconds  to 
microseconds  (ns  to  ~  20ps),  the  damage  mechanism  is  due  to  thermo-acoustic  effects 
[14],  In  this  region  the  diffusion  of  heat  is  negligible  and  the  damage  threshold  is 
fluence  dependent  (i.e.  J/cm2).  The  threshold  for  eye  damage  is  -0.1  J/cm2.  However, 
for  ultrashort  pulses  (subpicosecond),  the  data  shows  lower  threshold  damage  as  the  pulse 
width  decreases.  For  example,  the  damage  threshold  is  -0.01  J/cm2  for  picosecond  pulses. 
These  results  for  ultrashort  pulses  (-100  femtoseconds)  were  later  confirmed  by 
additional  experiments.  [18].  The  damage  mechanism  is  now  believed  to  be  due  to  self- 
focusing  and  laser  induced  breakdown  [19]. 

The  damage  threshold  [20]-[22]  also  depends  on  the  laser  spot  size  at  the  retina  where 
the  diffraction  limited  spot  size  is  22pm.  Measurements  of  the  damage  threshold  [16] 
were  performed  using  a  flashlamp-pumped  dye  laser  for  spot  sizes  ranging  from  22pm  to 
1400pm  (the  corresponding  field  of  view  is  from  1.5mrad  to  lOOmrad).  It  was  observed 
that  the  damage  threshold  increased  linearly  with  the  diameter  of  the  spot  size  in  this 
range.  However,  when  the  spot  size  exceeds  1400pm,  the  damage  threshold  is  fluence 
dependent,  i.e.  proportional  to  the  area  of  the  spot  [16].  Most  experiments  are  performed 
using  single  incident  pulses,  but  multi-pulse  experiments  have  also  been  performed. 
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Recent  experiments  [15]  have  investigated  multiple  exposures  with  repetition  rates  less 
than  several  kilohertz. 

Indirect  detectors  include  solid  detectors,  image  intensifiers  and  infrared  (IR)  night 
vision  detectors.  In  the  IR-II  and  IR-in  ranges  the  thermal  emissivity  of  various  objects 
is  used  to  produce  an  image  that  is  observed  by  an  IR  seeker.  Several  strategies  are  used 
to  damage  indirect  detectors.  For  example,  the  lens  in  front  of  a  typical  solid-state 
detector  will  focus  the  incoming  laser  beam  so  that  even  a  medium  energy  laser  can 
damage  the  detector.  Another  technique  for  sensor  damage  involves  destroying  the  front 
optics  of  the  detection  system. 

The  mechanism  responsible  for  optical  sensor  damage  also  depends  on  the  pulse 
duration:  below  Ins,  high  power  densities  can  be  reached,  leading  to  electrical 
breakdown;  above  Ins,  the  damage  originates  from  thermal  effects  that  are  fluence 
dependent  as  long  as  the  pulse  duration  is  shorter  than  the  thermal  time  constant  to 
become  power  density  dependent  [23], 

2  Mechanisms  of  Optical  Limiting 

Optical  limiting  was  first  observed  in  liquids  [12]-[13]  where  these  early  experiments 
demonstrated  that  the  limiting  could  be  due  to  self-focusing  or  self-defocusing  as  well  as 
absorption  and  scattering.  Later  experiments  showed  that  semiconductors  have  large 
optical  nonlinearities  [24].  Optical  power  limiting  experiments  conducted  in  CdS,  GaAs, 
and  CdSe  [25]  were  attributed  to  the  mechanism  of  two-photon  absorption.  Additionally, 
experiments  in  the  1960s,  showed  that  organic  dye  molecules  exhibited  optical  limiting 
property  [26],  This  was  recognized  as  reverse  saturable  absorption,  so-named  because  it 
is  due  to  the  large  absorption  cross-section  of  the  excited  state  of  the  molecule. 

During  the  last  several  years  many  mechanisms  used  for  passive  optical  limiting 
include  nonlinear  absorption,  nonlinear  refraction,  and  induced  scattering  [1]-[10].  The 
nonlinear  absorption  can  be  due  to  two-photon  absorption  or  reverse  saturable  absorption. 
Examples  of  nonlinear  refraction  are  the  electronic  Kerr  effect,  the  excitation  of  free 
carriers  or  optically  induced  heating.  Figure  1  shows  some  of  the  time  scales  for  the 
various  physical  phenomena. 

2. 1  Kerr  Nonlinear  Refraction 

At  high  input  energy  the  focal  point  of  the  beam  will  be  shifted  due  to  nonlinear 
refraction  [27],  therefore  the  energy  at  the  detector  will  be  limited.  Because  of  the  non- 
uniform  distribution  of  the  Gaussian  beam  profile,  the  nonlinear  medium  acts  as  either  a 
positive  or  a  negative  lens  depending  on  the  sign  of  the  nonlinear  refractive  index.  By 
proper  design  of  the  aperture  and  lens  structure,  the  nonlinear  medium  forms  a  lens  that 
can  increase  the  portion  of  the  beam  blocked  by  the  aperture. 

In  a  Kerr  nonlinear  medium,  the  change  of  index  is  due  to  3rd  order  nonlinear 
susceptibility. 

^Kerr  =  «27  0) 

For  a  Gaussian  beam,  the  transverse  intensity  distribution  is  given  by 
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(2) 


I  =  /0exp(-r2  lw2) 

The  refractive  index  distribution  is  given  by  [1] 

An  =  A n0  exp (-r2  /w2)  =  An0(l-r2  /w2) 


(3) 


An0=n2I0 

Another  critical  parameter  of  nonlinear  refractive  optical  limiter  is  the  position  of  the 
nonlinear  medium  in  the  device.  The  optimal  position  for  self-focusing  or  self-defocusing 
is  one  Rayleigh  range  before  or  after  the  focus  of  the  device,  respectively  which  can  be 
verified  by  Z-scan  experiments.  If  the  Kerr  nonlinear  medium  also  has  positive  dispersion 
then  pulse  splitting  occurs  as  the  pulse  propagates  [28]-[29],  These  effects  may  have  to 
be  taken  into  account  when  pulse  shaping  is  important. 


2.2  Nonlinear  Scattering 


Carbon-black  suspension  (CBS)  shows  broadband  nonlinear  scattering  effect  [30]- 
[31].  The  experimental  data  indicates  that  optical  limiting  due  to  CBS  is  fluence 
dependent,  and  effective  for  broadband  nanosecond  laser  pulses.  The  limiting  threshold 
of  CBS  is  around  100W  peak  power.  The  mechanism  of  CBS  optical  limiter  is  simple. 
The  microscopic-sized  carbon  particles  are  rapidly  heated  to  4000°  K  by  strong  linear 
absorption,  giving  rise  to  thermionic  emission,  which  in  the  presence  of  the  strong 
electric  field  leads  to  avalanche  ionization.  The  resulting  microplasma  then  rapidly 
expands  into  the  surrounding  liquid  and  strongly  scatters  the  incident  light  for  the 
duration  of  the  existence  of  the  plasma  (100ns).  In  addition,  since  the  carbon  particles  are 
ionized  during  the  process,  the  suspension  must  be  replenished  after  each  laser  exposure. 

2.3  Two  Photon  Absorption 


Two-photon  absorption  (TP A)  is  a  two-stage  process;  wherein,  the  electron  is  first 
excited  to  an  intermediate  state  (whose  lifetime  is  extremely  short),  followed  by 
absorption  of  a  second  photon  that  excites  it  to  the  final  state.  The  intensity  of  the  beam 
along  the  z-axis  can  be  obtained  by  solving  the  phase  and  amplitude  of  field  separately  in 
the  Maxwell  wave  equation.  If  the  excited  state  absorption  is  ignored,  then  the 
propagation  of  the  beam  intensity  is  given  by 

=  (4) 

oz 

where  a  is  the  linear  absorption  coefficient  and  /?2  is  the  TPA  coefficient  which  is 


related  to  the  imaginary  part  of  %(3),  the  third-order  nonlinear  susceptibility.  The 
solution  of  Eq.  (4)  is 


7(L)  = 


^ -  (5) 

(1+W) 

where  linear  absorption  is  assumed  to  be  zero.  It  is  easily  seen  that  the  light  transmission 
depends  on  the  initial  intensity,  TPA  coefficient,  and  thickness  of  the  medium.  An 
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increase  in  intensity  leads  to  a  decrease  in  transmission.  In  TPA  materials  it  is  useful  if 
the  first  excited  state  is  longer  than  the  pulse  width  so  that  each  electron  is  excited  once 
per  pulse.  However  a  short  upper  excited  state  is  required  to  reduce  the  saturation  of  the 
excited  state  absorption  that  diminishes  the  limiting  performance. 

2.4  Semiconductor  Nonlinearity 

Semiconductors  have  broad  absorption  bands  and  are  capable  of  producing  TPA  over  a 
large  wavelength  range  where  linear  absorption  is  low.  In  semiconductors  electrons  can 
be  excited  to  the  conduction  band  where  they  can  contribute  to  the  absorption.  In  the 
case  of  free  carrier  generation  in  semiconductors,  the  change  of  the  refractive  index  is 
given  by 

<6> 

2 nco  J-to 

where  neh  is  the  index  change  per  photo-generated  carrier  pair  [32]-[33],  However  a 

problem  with  semiconductors  is  that  they  tend  to  have  low  damage  thresholds  that  restrict 
their  applicability  in  optical  limiting.  The  typical  semiconductor  is  also  more  useful  for 
the  infrared  region  than  the  visible  region. 

2.5  Thermal  Nonlinearity 

The  effects  of  thermal  nonlinearity  depend  on  the  time  scale  involved.  The  refractive 
index  profile  of  the  medium  changes  as  a  light  pulse  passes  through  it  due  to  induced 
temperature  and  density  gradients.  This  mechanism  is  known  as  thermal  lensing  because 
the  refractive  index  change  follows  the  beam  shape  forming  a  lens-like  structure  in  the 
medium  [34],  For  steady-state  conditions  thermal  lensing  is  diffusion  dominated  [35], 
[36],  However  density  changes  resulting  from  the  temperature  gradient  are  described  by 
acoustic-waves  because  they  travel  at  the  speed  of  sound  [37]-[39]. 

When  the  laser  pulse  is  longer  than  a  few  microseconds  the  heating  effect  is  diffusion 
dominated  and  the  heating  of  the  material  is  described  by  the  diffusion  equation  [38] 

pc,~^2T  =  aLI  (7) 

ot 

where  p  is  the  density  of  the  medium,  Cp  is  the  specific  heat  at  constant  pressure,  k  is  the 
thermal  conductivity,  and  aLI  is  the  portion  of  the  laser  power  transformed  into  heat.  In 
nanosecond  region  the  effects  of  diffusion  can  be  ignored  and  the  induced  temperature 
change  for  any  spatial  position  is  given  by 

AT(t)  =  -—\aLI(t')dt'  (8) 

Pcp  0 

In  this  temporal  region,  the  thermal  effects  can  be  described  by  photon-acoustic  wave 
interaction 
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v2(a«)- 


1  32(A n) 


-V2  (AT) 


where  C]  is  the  velocity  of  sound.  This  expression  may  be  simplified  when  the  acoustic 
transit  time  is  less  than  the  pulse  width  and  is  given  by 

i  f  dtl  I  km  /i  /\\ 


where 


r*!  do 

UrJ.  2 n 

is  the  thermo-optic  constant.  However  for  picosecond  duration  pulses  both  the  thermal 
diffusion  and  photo-acoustic  effects  can  be  ignored. 

2.6  Reverse  Saturable  Absorption 

Reverse  saturable  absorption  (RSA)  occurs  when  the  excited  state  absorption  cross 
section  (cr^)  is  larger  than  the  ground  state  absorption  cross  section  (<jg ).  The  materials 

involved  are  usually  large  chromophores  [3]-[8],  The  process  is  modeled  by  several 
vibronically  broadened  electronic  energy  states.  In  general  a  five  level  system  is  involved 
as  shown  in  Fig.  2.  The  RSA  material  typically  consists  of  three  singlet  states  (ground 
state  and  two  excited  states)  coupled  with  two  triplet  states.  These  materials  have  a  large 
ratio  of  excited  state  to  ground  state  absorption  cross  section  and  can  achieve  a  large 
nonlinear  attenuation  while  maintaining  a  high  linear  transmittance  (>70%).  In  addition, 
RSA  materials  with  fast  singlet  excited  state  absorption  and  long-lived  triplet  state 
absorption  [40]-[43]  can  give  rise  to  optical  limiting  over  a  wide  range  of  pulse  widths 
(i.e.  nanoseconds  to  sub-picoseconds).  Modifications  to  the  molecular  structure  of  the 
chromophores  can  lead  to  enhanced  optical  limiting  in  specific  spectral  ranges. 
Consequently  theoretical  [44]-[48]  and  experimental  [3]  studies  are  in  progress  in  several 
laboratories. 

In  the  simplified  case  [1],  [7]  the  transmission  through  the  RSA  materials  is 
described  by  the  fluence  (F)  equation  given  by 

ey^_  F(0)exp(-aL) 


F(L)  = 


1  +  {gjx.  /  2ha))F(0)L 


=(l-exp(-aZ))/a 

where  F(0)  is  the  incident  fluence.  Similar  to  the  TP  A  process,  the  transmittance  of  RSA 
depends  on  the  incident  fluence,  excited  states  absorption  cross-section,  linear  absorption 
coefficient  and  effective  thickness. 


3.  RSA  Materials 
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A  number  of  materials  are  used  for  optical  limiting  such  as  fullerenes,  porphyrins, 
phthalocyanine,  naphthalocyanines  and  neutral  nickel  dithiolenes  [6],  Some  of  the  main 
optical  limiting  candidates  are  described  next. 

3 . 1  Fullerenes-Carbon  60  /  Carbon  70 


Carbon  60  (C60)  has  excited-state  absorption  greater  than  the  ground  state  from  -400 
to  900nm,  which  covers  the  entire  visible  spectrum.  Therefore  C60  is  a  broadband 
optical  limiting  material  [49]-[57],  Due  to  a  long  triplet  excited-state  lifetime  >40ps,  C60 
can  protect  against  pulse  duration  of  hundreds  of  microseconds  pulses.  The  performance 
of  C60  optical  limiter  also  depends  on  host  materials  [49],  e.g.  toluene  solution, 
polymethyl  methacrylate  (PMMA)  and  solid  xerogel  matrices. 

Carbon  70  (C70)  has  a  smaller  ratio  of  excited-state  to  ground  state  absorption  and  a 
higher  threshold  for  optical  limiting  [52],  [53].  In  the  visible  wavelength  range,  C70  has 
stronger  absorption  than  C60  because  of  the  reduced  symmetry  of  C70  relative  to  C60 
and  a  lower  triplet  state  quantum  yield.  Therefore  C60  and  C70  have  different  clamping 
intensities  and  the  limited  output  intensities.  The  values  of  the  mixture  C60/C70  solution 
are  higher  than  that  of  C60  solutions  and  lower  than  that  of  C70  solutions.  Variable 
limited  output  intensity  can  be  obtained  by  changing  the  C70  content  in  the  solution  of 
C60/C70  mixture  or  the  concentration  of  the  solution  [51]. 

3.2  Porphyrin  and  related  derivatives 

Optical  limiting  in  porphyrin-based  materials  was  first  reported  in  1985  [6], 
Thereafter  an  extensive  range  of  porphyrins  and  their  expanded  family  have  been 
intensively  studied  [58]-[73].  Figure  4  shows  a  basic  porphyrin  structure.  To  date  the 
largest  ratio  of  excited-state  to  ground-state  cross  section  absorption  reported  is  about  a 
factor  of  45  at  532nm  [58]  in  tetra  (trimethylsilyethynyl)  porphyrin  (TTMSAP(Pb)). 
Preliminary  work  shows  porphyrins  with  closed-shell  metal  atoms  have  relatively  long 
excited  state  lifetimes.  A  recent  halogenation  study  shows  that  bromine  atoms  substituted 
in  porphyrins  [59]  increase  the  triplet  yield  via  the  heavy  atom  effect  and  decrease  the 
triplet  lifetime.  One  molecular  design  strategy  is  to  maintain  a  planar  and  symmetric 
structure  to  ensure  that  the  spectra  of  linear  absorption  features  are  narrow.  Ab  initio 
electronic  structure  calculations  are  used  to  examine  the  structures,  excited  state  spectra, 
and  ionization  potentials  of  porphyrins  [6],  [44]-[48]. 

Recent  experimental  work  has  shown  that  heavy  central  atoms,  halogens,  and  electron 
donating  groups  in  porphyrins  can  enhance  RSA  and  optical  limiting.  However,  the 
mechanisms  involved  in  enhancing  optical  limiting  are  not  understood.  Nguyen  et  al.[60] 
studied  the  effects  of  zinc  and  halogens  in  porphyrins,  in  order  to  examine  their  possible 
role  in  enhancing  RSA,  using  ab  initio  electronic  structure  calculation  and  configuration 
interaction  method  using  single  substitutions  (CIS)  calculation.  The  CIS  calculations 
examined  the  effect  of  the  central  metal  and  (3  -halogenation  on  the  electronic  spectra  of 
porphyrins.  Their  calculation  of  excitation  energy  agrees  well  with  experiments.  Thus 
these  computational  methods  may  provide  an  approach  for  modeling  the  properties  of 
optical  limiters  on  a  molecular  scale  [60], 
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4  General  five  level  model  for  RSA 
4. 1  Rate  equations 

The  well-known  phenomenological  rate  equations  [6]  describing  the  electronic  state 


populations  for  the  five-state  model  are 

8NJ&  =  -cJNJho),  +kl0N{  +k30N3  (13) 

8N |  /  dt  a0lI  N, o  /  ho)  q  (ct X2J  I  ho)  ^  ^io  ^13)^1  ^21^2  (14) 

dN2/dt  =  an!NJho)0-k2XN2  (15) 

dN3  /  dt  —  —(cr34I  /  ho)  0  +  k3Q)N j  +  kX3Nx  +  k43N4  (1  6) 

dN4  ldt  =  a34  I N3 /h(D0  -  k43N4  n  7x 


where  Np  is  the  electron  number  density  of  the  state  p,  apq  is  the  absorption  cross-section 
for  electron  pumping  from  the  state  p  to  the  state  q  and  kpq  is  the  decay  rate  from  the  state 
p  to  the  state  q.  The  assignment  of  the  electron  densities,  Np,  in  Fig.  2(e)  is  as  follows: 
No  corresponds  to  the  ground  level  of  the  singlet  state.  So;  Ni  corresponds  to  the  first 
excited  level  of  the  singlet  state,  Si;  N2  corresponds  to  the  second  excited  level  of  the 
singlet  state,  S2;  N3  corresponds  to  the  first  excited  level  of  the  triplet  state,  Ti;  and  N4 
corresponds  to  the  second  excited  level  of  the  triplet  state,  T2.  In  these  materials,  the 
total  number,  Nt,  of  electrons  is  conserved  such  that  No  +  Ni  +  N2  +  N3  +  N4  =  Nt.  The 
initial  population  is  given  by  No=Nj  and  Ni=N2=N3=:N4=0. 

4.2  Intensity  Equations 

The  optical  pulse  has  an  intensity  I  with  angular  frequency  ©0,  and  h  is  Planck’s 
constant.  Using  a  density  matrix  approach,  we  [74]-[77]  obtained  the  induced  electric 
polarization  vector  from  the  dipole  moment  perturbation  calculation  for  the  RSA  material 
and  reduced  the  general  propagation  equation  in  our  model  to  an  equation  for  the  pulse 
intensity  that  is  in  agreement  with  previous  results.  This  reduced  propagation  equation 
for  the  pulse  intensity  is  given  by 

^  =  -[(T0,  AT,  +  <7„  AT,  +  CTmAT4  ]/  (18) 

To  date  most  studies  of  RSA  have  been  experimental,  with  a  few  theoretical  studies 
using  the  reduced  intensity  model.  More  complete  theoretical  studies,  including 
propagation  effects,  are  needed  in  order  to  guide,  validate  and  reduce  the  number  of 
experiments. 


8 


For  most  applications,  the  RSA  material  is  contained  in  a  matrix  substrate.  In  certain 
cases,  the  matrix  substrate  can  exhibit  Kerr  nonlinearity  and  material  dispersion.  The 
polarization  vector  for  this  material  can  be  expressed  as 

P  =  PL  +  P„L.  09) 

where  Pl  is  the  linear  polarization  vector  and  Pnl  is  the  nonlinear  polarization  vector. 
The  nonlinear  polarization  term  can  be  written  as 

Pnl=sosnlE »  (2°) 

where  E  is  the  electromagnetic  field  and  Snl  is  the  nonlinear  contribution  to  the  dielectric 
constant.  It  is  described  by 

sNL=\x('l\Ef,  (21) 

where  ^3).  is  the  third  order  susceptibility.  The  index  of  refraction  contains  both  a  linear 
frequency  dependent  term  and  a  nonlinear  intensity  dependent  term  defined  by 

n(a,\E\2)  =  n(to)+n2\E\\ 


(22) 


The  electromagnetic  field  E(r,z,t)  propagating  in  a  matrix  substrate  containing  the 
RSA  material  is  written  in  terms  of  a  slowly  varying  envelope  A(r,z,t)  and  the  carrier 
frequency  oo  of  the  electromagnetic  field,  such  that 

E(r,z,t)  =  A(r,z,t)e1k'*-ia,°t  +c.c.  (23) 

Using  the  slowly  varying  envelope  approximation,  we  obtain  the  equation  for 
propagation  of  the  electromagnetic  field  in  the  matrix  substrate  containing  the  RSA 
material  [74]-[75] 


{l+*'l]42>r'')+5;vi'4(2’r'')4iji'4z’r'')+[5fi'1'415. 

+^mN0(z,t]Af)+GnNx(zj\A^  +  GMN^z,t]A^A{z,r,t)  =  0 


A{z,r,t) 


(24) 


where  V]_ 


Id  ,  8nk 
S  +  Z~’  K=  — 


and  the  last  term  describes  the  RSA  material 


dr  r  dr  dm" 

whose  intensity  dependence  is  given  in  Eqs.  (13-17).  Equation  (18),  in  the  absence  of 
diffraction,  dispersion  and  Kerr  nonlinearity,  reduces  to  the  intensity  equation.  It  should 
be  noted  that  the  Kerr  nonlinearity  is  real  and  conserves  energy.  In  the  absence  of  RSA, 
but  together  with  anomalous  dispersion,  the  Kerr  nonlinearity  can  give  rise  to  localized 
solitary  waves  and  solitons.  However,  the  RSA  nonlinearity  is  pure  imaginary  and  leads 
to  an  intensity  dependent  energy  loss.  Because  of  their  nonlinear  dependence  on 
intensity,  RSA  materials  behave  very  differently  than  linear  absorbers.  Their  ideal 
feature  is  the  ability  to  transmit  low  intensity  light,  but  absorb  high  intensity  light. 


9 


We  assume  an  incident  pulse  of  the  form,  A(r,0, t)  =  exP 


r-r2' 

K2al. 


exp 


( 


,  where 


ao  is  the  spatial  beam  waist,  To  is  the  1/e  half-width  of  the  temporal  intensity,  and  Ao2  is 
the  incident  peak  intensity.  We  use  the  transformation  [75], 


n-—  n-—  P--J—  r 

Q  k<flyr  t,  ' 


V  a„ 

to  obtain  the  non-dimensioned  equation 


(25) 


.  a  l-j  i  a2 
»«+t  v2--r- 


4  '  2  dr2 


Q{£,p,  *)  +qp|2(£  a  *)|20(£>a  7) 


(26) 


^<70,JV0(ft  r,|fif )+ <712Ar,(*  r.|ef )  +  crMiVs(f,  r,|g|2)]e(f,  p,r)  =  0 
coupled  to  the  rate  equations 


i 

o 

K 

0 

*30 

0 

X" 

Ai 

-(A2+* io+*i3) 

*21 

0 

0 

JV, 

n2 

= 

0 

A 

A2 

*21 

0 

0 

^2 

0 

^13 

0 

o 

+ 

i 

*43 

W 

0 

0 

0-34 

—  *43  j 

IxJ 

(27) 


where  ap  =  (i^J  ,r  - =  <r„  4 |gf  /ha„  and 

W  denotes  differentiation  with  respect  to  time.  Equations  (26)  and  (27)  represent  the  set 
of  coupled  equations  that  must  be  solved  numerically. 

The  numerical  procedure  combines  a  split  step  beam  propagation  method  for  the 
dispersive  terms  with  a  finite-difference  method  for  the  diffractive  terms.  The  nonlinear 
terms  are  treated  in  the  time  domain.  The  numerical  solution  is  solved  in  the  form 


10 


|2  =  (/>„+/>„+«!)& 


D  =  —V2 

■^DF  “  4  V’ 

NL  =  NLl+NL2 
NL\  =iap\Qf 

F4=-^0,F„«,r,|er)+<T,!F,«,r,|(3|2)+Cr!4Ar!(f,r,|e|2)] 


(28) 


where  Dos  (Ddf)  is  the  dispersion  (diffraction)  operator,  NLi  describes  the  Kerr 
nonlinearity  operator  and  IsOU  describes  the  nonlinear  absorption  operator  resulting  from 
the  RSA  material.  For  the  symmetrical  split  beam  method  (BPM),  we  write  the  equation 
as 


(29) 


where  the  integral  is  approximated  by 


AGE,  (f  ^  [JVZ,  (f + A£)  +  M,  (ft ] . 


(30) 


5.  Numerical  Method 

A  schematic  diagram  of  the  numerical  code  processes  is  given  in  Fig.  3.  Some  of  the 
material  parameters  include  the  cross-sections  for  absorption  (cr0l,crn,<J34 ),  the  decay 

rates  for  the  five-level  model  ( ,  k4% ,  kll ),  the  carrier  concentration,  the  Ken- 

nonlinear  coefficient,  the  2-photon  absorption  coefficient,  the  group  velocity  dispersion, 
the  linear  index  of  refraction  and  linear  absorption.  Many  of  these  parameters  can  be 
varied  in  the  numerical  calculations  to  optimize  the  performance  and  compared  with  ab 
initio  quantum  molecular  calculations  to  guide  molecular  modeling.  Examples  of  input 
optical  beam  parameters  include  pulse  intensity,  width,  radius,  shape,  and  repetition  rate. 
Analysis  of  the  output  files  can  identify  any  changes  in  pulse  reshaping  that  may  affect 
nonlinear  transmission.  In  other  words,  pulse  reshaping  such  as  self-focusing  or  de- 
focusing  or  temporal  compression  or  expansion  may  result  in  beneficial  optical  limiting 
results.  Some  device  parameters  include  the  number  of  lens  and  their  focal  length  and 
distance  from  each  object  and  the  number  of  elements  in  the  optical  limiter.  The 
numerical  calculations  can  be  used  to  optimize  the  device  design  and  thus  minimize  the 
number  of  costly  experiments  that  have  to  be  performed. 
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The  numerical  procedure  combines  the  split  step  beam  propagation  method  for  the 
dispersive  terms  and  the  finite-difference  method  for  the  diffractive  terms.  The  nonlinear 
terms  are  treated  in  the  time  domain.  Numerical  calculations  are  performed  with  non- 
dimensioned  equations. 

5.1  Numerical  Method  and  Operator  Equations 


The  mathematical  terms  [74]-[76]  are  evaluated  using  a  discrete  space  and  time  grid 
given  by  r  -*  r(.  p  -»  p}  and  the  field  can  be  written  in  the  shorthand  notation, 

Q"j.  The  dispersion  and  diffraction  terms  are  calculated  in  the  frequency  domain  using 
the  fast  Fourier  transform  method;  however,  the  nonlinear  terms  are  calculated  in  the  time 
domain. 


In  the  frequency  domain,  the  dispersion  operator  becomes 


DDS(p)  = 


lyco 2 


The  linear  diffraction  term  is  given  by 

^  =  kV2pQ{Z,p,co),  k='~. 
This  term  is  solved  by  the  finite-difference  scheme, 

s_q  ,gy-g, 

a#  A<f  ’ 
and  the  Crank-Nicholson  method 

ne= -erji +e:,,, -sv.3 


+0% ~2&j  +0u-i] 


Ap  =  pj+l-pj 

for  0  <Pj  <  Pmax- 

Combining  Eqs.  (33)  and  (34),  the  field  at  /  2  is  given  by 


(31) 


(32) 


(33) 


(34) 


k  A  £ 


Ap  2 


1 


1 


2  pj  Ap 


\QijU  + 


tcAi; 

(ApF 


+1 


\Qu 


K  A£ 


Ap  2 


(  1  l' 
+ 


2  Pj  Ap 


On+l  = 


(35) 


K  Ai f 
~~Ap~2 

However  at  p=0, 


2  Pj  Ap 


\Q"h  ~ 


'  *  At; 


(AP)2 


CC/  + 


k  A£ 


) 


Ap  2 


1  1 

-+- 


2  pj  Ap 


Ql 


j+\ 
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(1 + <p)QSl  -<mx  -  (1  ■ -  <p)Qlo  -  <pQl  =  o, 

2k  A£ 

<P  = - 

P\ 

At  the  maximum  radial  value,  (p  =  /?max  =  pN)  we  obtain 

e  a".,  +  (i  -  eyns + 90V-,  -  ( i+ ejec. = o. 


(/-V  P.V-1 ) 

Using  Eqs.  (35)-(37),  Eq.  (33)  is  written  in  the  form 


«•=(&. . 

and  Ei  and  E2  are  tri-diagonal  matrices.  The  matrix  Ei  is  composed  of  the  sub-diagonal, 
diagonal  and  super-diagonal  elements, 


kA{ 

1 

$ 

+ 

1 

£(*>,)  =  . *^1- — - — +— - — L.,1-0, 

Pj+\  Pj~\  [Pj+\  Pj  Pj  Pj- ij 


\  AT  1  1 

CK)  =  “?>•••» - - — i  — + - -N-»° 

_  Pj+x-Pj-x  [2Pj  Pj+X~Pj  J  J 

and  the  matrix  E2  is  composed  of  the  sub-diagonal,  diagonal  and  super-diagonal  terms 
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<*K)= 


o,...  - 


*rA<f  1  1 


P;+1-P;-.  [2py  Pj-p,-, 


,-0 


*K)  = 


irA£ 


P./+1  P;-l 


1 


1 


p j+1  pj  p  j  p j-\  J 


,...,1+0 


(40) 


/K)= 


tcAg 


PJ+\  Pj~\ 


1  1 


-+- 


P,+i"P,J 


■,...,o 


5.2  Propagation  Equation 


The  propagation  including  the  RSA  nonlinearity  is  more  complex  because  of  the 
coupling  between  the  propagation  equation  and  the  rate  equations.  The  propagation  term 


for  the  RSA  is  written  as  the  coupled  set  of  equations 

^  =  -5(o-0l^+<71!#,+<TJ4ArJ)e,  (41) 

An(1  t,  |ef  )  =  M[l  t,$)n((,  r,|gf),  (42) 

The  propagation  can  be  expressed  as 

er'=exp[-5^(off)]a",  (43) 

(jN  =  ^"oi  1/2,  ^1-1/2)  "^0 (^«+l/2, ^"i+l/2  )] 

+  ^12  [^l(^«+l/2,  ^j-1/2  )  +  (^w+1/2,  ^i+l/2  )]  (^) 


»+l/2,7M/2  )  +  -^3  (^»+l/2, 7i+l/2  )] 


However,  the  N’s  must  be  determined  from  the  rate  equation,  and  depend  on  the  intensity 
|(2[2  •  The  rate  equations  can  be  solved  in  matrix  form  as 


fl+1/2  ?  7j+1/2 


A  T 


m\[n& 


n+ 1/2 »  7i+l/2 


1 )  +  N{£„+\/2  >  ^1-1/2)] 


(45) 


Collecting  like  terms  gives  the  matrix  expression 

i-i 


N(€„+ 1/2, 7,+1/2)  — 


At  a 

1-—M 

2 


At  * 
1  +—M 
2 


-^(^1+1/2  >  )  > 


(46) 


expanding  the  inverse  matrix  and  keeping  terms  to  second  order, 
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(47) 


<V(C„! ,rM „ )  =  [l  +  A  r  M + ( Ar  /  2  »(£,„„ ,  r, 

Equation  (47)  is  used  in  the  propagation  equation. 

6.  Results  and  Conclusions 

The  following  sections  highlight  the  main  results  and  conclusions. 

6.1  Initial  Results 

Some  initial  results  are  described  in  this  section.  Figure  4  shows  the  pulse  energy  as  a 
function  of  space  and  time  for  an  incident  energy  of  10  pJ  at  four  propagation  distances 
(z=0.28  cm,  0.57  cm,  0.86  cm,  and  1.15  cm).  The  pulse  shapes  are  plotted  as  functions  of 
time  and  radius,  and  the  contour  plot  of  the  pulse  is  given  at  the  top  of  each  pulse  shape. 
The  negative  (positive)  times  correspond  to  the  leading  (trailing)  edge  of  the  pulse.  It  can 
be  seen  from  Fig.  4,  that  the  trailing  edge  of  the  pulse  becomes  less  steep  as  the  pulse 
propagates  and  the  peak  of  the  pulse  advances  in  time  relative  to  the  moving  coordinate 
frame.  For  plotting  purposes  the  pulses  are  normalized  to  a  height  of  one  unit.  For 
convenience  the  values  of  the  peak  energy  at  each  propagation  distance  is  given  in  Fig.  4. 
These  values  show  the  decrease  in  energy  as  the  pulse  propagates.  While  the  overall 
pulse  shape  exhibits  small  changes  as  it  propagates,  the  energy  decreases  rapidly.  An 
example  of  the  carrier  density  is  shown  in  Fig.  5.  For  convenience,  the  carrier  densities 
are  normalized  by  the  expression 

n0  +«i  +n2  +«s  +n4  =  1.0, 

nt  =  NJ  Nt  V  ' 

where  l  =0,1, 2, 3, 4  and  the  initial  condition  of  the  carrier  density  before  the  incidence  of 
the  pulse  is  no=l  .0  and  m=n2=n3=n4=0.  In  all  cases  it  was  found  that  the  value  of  iu  was 
negligible;  therefore,  it  was  not  included  in  the  plots  with  the  other  carrier  densities. 
Figure  5  shows  plots  of  the  carrier  densities  at  z=0.28  cm.  The  maximum  values  of  the 
carrier  densities  are  given  in  the  figure  caption  as  well  as  the  minimum  value  for  no.  The 
minimum  value  for  all  other  n’s  is  zero,  so  it  is  not  necessary  to  indicate  their  values. 
Initially  the  entire  carrier  population  is  in  So  (i.e.  no  =1.0),  but  as  the  leading  edge  of  the 
pulse  enters  the  RSA  material,  carriers  are  excited  from  the  ground  state  (So)  to  the  first 
excited  state  (Si).  Because  the  entire  carrier  density  is  originally  in  So,  the  pulse 
absorption  is  determined  by  the  product  <70,/70 1  at  the  leading  edge  of  the  pulse.  As  the 
first  excited  state  becomes  populated  (i.e.  ni>0),  the  absorption  of  the  pulse  is  determined 
by  the  relationship  (&01n0  +<yunx)I.  The  population  of  the  second  excited  state  (i.e.  n2) 
decays  very  rapidly  because  k2f 1  <  r0 .  The  decay  rates  kio'1  and  kn'1  >  r0  so  the  first 
excited  state  (i.e.  Si)  can  maintain  some  population  (i.e.  ni>0)  for  a  time  greater  than  r0 . 
However,  as  shown  in  Fig.  5,  there  is  minimal  decay  to  the  triplet  state  Ti.  Therefore, 
even  though  <734  is  about  the  same  size  as  012,  there  is  little  population  of  T2  (i.e.  n4~0) 
which  is  verified  by  the  numerical  calculations.  Because  Si  remains  populated  for  a  time 
greater  than  r0,  the  trailing  edge  of  the  pulse  intensity  experiences  greater  absorption 
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than  the  leading  edge,  thus  causing  the  peak  of  the  pulse  to  move  forward  in  time  relative 
to  the  moving  frame. 

Figure  6  shows  the  pulse  energy  as  a  function  of  space  and  time  for  an  incident  energy 
of  500  pj.  In  this  case,  new  details  of  the  pulse  propagation  behavior  are  observed  using 
the  space  and  time  dependent  code.  Initially  strong  carrier  depletion  occurs  at  the  center 
of  the  pulse  resulting  in  a  slower  decay  for  the  peak  energy.  This  behavior  is  in  contrast 
with  the  previous  results  shown  in  Fig.  4.  In  order  to  investigate  the  origin  of  this 
behavior  the  pulse  shapes  are  plotted  at  four  propagation  distances  (the  same  distances 
used  in  Fig.  4).  Figure  6  shows  a  new  lobe  developing  on  the  leading  edge  of  the  pulse 
(z=0.57  cm).  As  the  pulse  propagates  the  high  intensity  region  is  absorbed  by  the 
material  leading  to  pulse  splitting  at  z=  1.1 5  cm.  Figure  7  shows  the  carrier  densities  at 
z=0.28  cm.  The  leading  edge  of  the  pulse  excites  the  carriers  from  the  ground  state  (So) 
to  the  first  excited  state  (Si).  However  because  of  the  high  pulse  energy.  So  is  almost 
completely  depleted  (i.e.  no  decreases  to  a  minimum  value  of  ~10'4=  0).  At  the  leading 
edge  of  the  pulse,  the  carrier  density  is  increased  to  the  first  excited  state  Si  (i.e.  ni>0). 
However  at  the  center  of  the  pulse  the  energy  is  sufficient  to  populate  the  second  excited 
state  S2  (i.e.  n2  reaches  a  maximum  value  of  ~0.8).  This  maximum  population  of  S2 
occurs  only  near  the  pulse  center  because  the  decay  rate  for  this  state  is  faster  than  the 
temporal  pulse  width.  In  fact,  this  rapid  population  and  depopulation  of  S2  actually 
results  in  a  depletion  of  ni  near  the  pulse  center  as  shown  by  the  contour  plot  of  ni  in  Fig. 
7.  However,  the  trailing  edge  of  the  pulse  still  has  sufficient  energy  to  create  significant 
population  of  Si  (i.e.  m  reaches  a  maximum  value  of  0.987  for  temporal  values 
corresponding  to  the  trailing  edge  of  the  pulse)  as  shown  in  Fig.  7.  The  decay  rate  from 
Si  to  So  is  longer  than  the  temporal  pulse  width  so  Si  does  not  decay  rapidly.  The  state 
Ti  is  not  heavily  populated  because  of  the  relatively  long  decay  rate  of  Si  to  Ti.  (ty 
reaches  a  maximum  value  of  ~10‘2)  Numerical  calculations  show  that  114  =  0.  As  the 
pulse  propagates  (z=0.57  cm)  through  the  RSA  material  a  temporal  lobe  forms  on  the 
leading  edge  of  the  pulse  as  shown  in  Fig.  6.  This  temporal  lobe  is  caused  by  the  complex 
interplay  between  the  pulse  intensity,  absorption  cross-sections  and  decay  rates  of  Si  and 
S2.  As  the  pulse  propagates  to  z=0.86  cm,  the  higher  intensity  portion  of  the  pulse  is 
reduced  relative  to  the  lobe  on  the  leading  edge  because  of  stronger  absorption.  The 
carrier  density  remains  mostly  in  So  because  the  intensity  at  this  distance  is  only  0.018 
which  is  not  sufficient  enough  to  cause  much  carrier  excitation.  Then  as  the  pulse 
propagates  to  z=1.15  cm  temporal  splitting  develops  because  the  trailing  edge  of  the 
pulse  has  been  more  strongly  absorbed  owing  to  its  greater  intensity. 

6.2  Three  Distinctive  Regions 

Numerically  we  observed  three  distinct  regions:  (1)  linear  absorption  due  to  carrier 
excitation  from  the  ground  state  to  the  first  excited  state;  (2)  limiting  region  due  to  carrier 
excitation  from  the  first  excited  state  to  the  second  excited  state;  and  (3)  saturation  region 
at  high  incident  energy.  In  the  linear  region  (low  incident  energy  <  10  uJ)  the  absorption 
is  due  to  the  ground  state  to  the  first  excited  singlet  state  and  there  is  a  symmetric 
reduction  in  the  pulse  shape.  This  region  is  not  shown  here  due  to  limitations  of  space. 
Figure  8  shows  the  2-D  and  3-D  intensity  and  carrier  density  plots  in  the  limiting  region 
(50  uJ).  Part  (I)  of  this  figure  shows  the  2-D  plots  of  the  pulse  Intensity  and  the  carrier 
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density;  and  part  (II)  shows  the  3-D  plots  of  the  carrier  densities.  In  this  case,  the 
absorption  is  due  to  the  transition  from  the  first  excited  singlet  state  to  the  second  excited 
singlet  state  leading  to  greater  absorption.  The  leading  edge  of  the  pulse  excites  the 
carriers  from  the  ground  state  to  the  first  excited  state  and  the  trailing  edge  of  the  pulse  is 
absorbed  by  the  carrier  excitation  from  the  first  excited  state  to  the  second  excited  state. 
As  a  result,  most  absorption  occurs  in  the  trailing  edge  of  the  pulse  resulting  in  a 
temporal  shift  of  the  pulse  toward  the  leading  edge. 

Figure  9  shows  the  2-D  and  3-D  intensity  and  carrier  density  plots  in  the  saturation 
region  (1000  uJ).  Part  (I)  of  this  figure  shows  the  2-D  plots  of  the  pulse  intensity  and  the 
carrier  density;  and  part  (II)  shows  the  3-D  plots  of  the  carrier  densities.  At  high  input 
energy  (-1000  uJ)  the  carriers  are  excited  to  the  second  excited  state  at  the  leading  edge 
of  the  pulse  resulting  in  less  absorption  during  the  pulse  duration;  thereby,  increasing  the 
transmission  as  the  input  energy  increases. 

Further  numerical  results  show  that  variation  of  the  material  parameters  such  as  the 
decay  rates,  the  cross  section  of  absorption  of  the  material,  and  initial  carrier  density,  in 
addition  to  the  pulse  properties  such  as  the  pulse  width  and  spatial  width  can  enhance  the 
limiter  performance.  Therefore  a  detailed  microscopic  understanding  of  the  pulse  and  the 
carrier  densities  in  each  of  the  three  physical  regions  is  required  to  describe  the  dynamics 
of  limiter  behavior. 

6.3  Propagation  Through  RSA  Material  in  Kerr  Medium 

In  general  the  RSA  material  is  embedded  in  a  matrix  substrate,  which  may  exhibit  the 
Kerr  nonlinearity.  This  nonlinearity  may  lead  to  self-focusing  of  the  light  beam  as  it 
propagates  through  the  material.  Figure  10  shows  a  plot  of  the  transmission  as  a  function 
of  the  coefficient  of  the  Kerr  nonlinearity  and  demonstrates  that  self-focusing  can 
enhance  the  nonlinear  absorption  of  the  RSA  material.  Because  the  linear  absorption  is 
quite  significant  in  the  RSA  material  much  of  the  laser  energy  is  absorbed  before 
significant  self-focusing  can  occur.  Therefore  we  investigate  the  effects  of  a  two  layer 
absorber  in  which  the  first  layer  exhibits  the  Kerr  nonlinearity  and  the  second  layer 
exhibits  RSA  absorption.  First  we  describe  various  pulse  shaping  techniques. 

6.4  Pulse  shaping  techniques 

Examples  of  pulse  shaping  with  nonlinearity  are  shown  in  Figs.  11,  12  and  13. 
Pulse  shaping  techniques  include  reshaping  the  pulse  in  the  time  domain  i.e.  pulse 
compression  or  pulse  broadening  and  reshaping  in  the  spatial  domain  i.e.  self-focusing  of 
self-defocusing. 

6.5  Spatial  pulse  shaping:  Self-focusing 

It  has  been  noted  previously  that  for  RSA  limiters,  the  material  is  placed  in  a  focal 
plane.  However  this  can  be  cumbersome  for  some  applications.  Another  approach 
involves  using  a  Kerr  nonlinear  medium  before  the  RSA  material  to  self-focus  the  light. 
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Figure  16  shows  the  use  of  a  Kerr  nonlinear  medium  to  self-focus  light  into  a  RSA 
tandem  limiter.  Figure  1 1  shows  change  in  transmission  in  an  optical  limiter  when  a  Kerr 
nonlinear  medium  is  used  before  the  nonlinear  absorber  material.  This  figure  shows  the 
transmission  as  a  function  of  incident  energy. 

6.6  Temporal  pulse  shaping  techniques 

In  the  time  domain,  pulse  shaping  can  be  accomplished  by  using  Kerr  nonlinearity 
and  GVD  (negative  GVD  leads  to  pulse  compression  and  positive  GVD  leads  to  pulse 
expansion).  For  example,  pulse  shaping  may  be  used  before  the  RSA  material  in  order 
to  change  nanosecond/microsecond  pulses  to  picosecond  pulses  or  vise  versa.  This  can 
then  be  used  to  take  advantage  of  the  decay  rates  and  kinetics  of  the  RSA  material. 
Additionally  certain  detectors  may  be  more  sensitive  to  broad  temporal  pulses  and 
some  may  be  more  sensitive  to  short  temporal  pulses.  Thus  pulse  reshaping  using 
passive  nonlinearity  in  a  tandem  optical  limiter  may  provide  a  significant  advantage 
that  would  not  be  possible  without  the  Kerr  nonlinear/dispersive  medium. 

6.7  Self-focusing  and  RSA 

The  effects  of  the  combination  of  self-focusing  and  RSA  are  shown  in  Figs.  14- 
16.  The  input  parameters  are  To=  4.0  ps,  Ro  =45um,  Lambda=  532  nm  and  carrier 
density=  6.0322e-4  nm'3.  Figure  14  shows  the  pulse  compression  after  one-diftffaction 
length  in  a  Kerr  nonlinear  medium  where  the  peak  intensity  increases  by  a  factor  of  18. 
Figure  15  shows  the  propagation  of  the  compressed  optical  pulse  in  0.7  mm  of  RSA 
material.  It  shows  the  depletion  of  the  ground  singlet  state  and  the  population  of  the 
second  excited  singlet  state.  There  is  veiy  little  population  of  the  triplet  state.  Figure  16 
shows  an  increase  in  absorption  when  self-focusing  is  used  to  precede  the  RSA  material. 
This  material  acts  as  a  nonlinear  lens. 

6.8  Multiple  Pulses 

Substantial  efforts  have  been  used  to  investigate  transmission  for  various  optical 
limiters  using  a  single  incident  pulse.  However,  lasers  often  produce  a  train  of  pulses. 
We  investigated  numerically  the  effects  of  multiple  pulses  and  pulse  separation.  The 
input  condition  for  multiple  pulses  is  given  by 

Q(0,r)  ^^exP~\(T~krr)2 /2to]  (49) 

k=Q  L  J 

where  rris  the  pulse  separation  and£0  is  the  pulse  energy.  Figure  17  shows  a  schematic 
diagram  of  the  multiple  pulses.  Figure  18  shows  the  transmission  curve  for  a  train  of 
seven  pulses  as  a  function  of  input  energy  for  pulse  separations  of  xr  =2  and  10.  Our 
results  show  that  for  closely  spaced  pulses  (xr  =2)  there  is  a  reduction  in  transmission  in 
the  region  of  excited-state  absorption. 

Detailed  analysis  of  the  pulse  shapes  indicate  that  this  effect  is  due  to  pulse 
coalescence  that  results  in  increased  input  energy.  Figure  19  shows  2-D  and  3-D  plots  of 
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the  pulse  train  intensity  for  seven  incident  pulses.  Fig.  (a)  and  (b)  show  the  pulse  train 
before  the  excited-state  material.  Fig.  19(b)  shows  that  the  seven  incident  pulses  have 
merged  into  five  pulses  and  their  peak  energy  has  increased  by  about  60%.  Fig.  19  (c) 
and  (d)  show  the  pulse  train  after  transmission  through  the  excited-state  material.  The 
first  pulse  in  the  train  excites  carriers  to  the  excited  state  absorption  level;  however,  the 
carriers  do  not  have  sufficient  time  to  return  to  the  ground  state  before  the  next  pulse 
arrives.  Eventually,  an  equilibrium  state  among  the  carrier  densities  is  achieved  such  that 
the  trailing  pulses  reach  constant  transmission  energy. 

We  also  investigated  the  pulse  shapes  for  the  case  of  large  pulse  separation  (xr  =  10). 
Figure  20  shows  2-D  and  3-D  plots  of  the  pulse  train  intensity  for  seven  incident  pulses. 
Fig.  (a)  and  (b)  show  the  pulse  train  before  the  excited-state  material.  Fig.  20(b)  shows 
that  the  seven  incident  pulses  do  not  merge  and  their  peak  energy  does  not  increase. 
Therefore  Fig.  1  shows  greater  transmission  for  this  input  energy  than  for  the  closely 
spaced  pulse  train.  Fig.  20(c)  and  (d)  show  the  pulse  train  after  transmission  through  the 
excited-state  material.  Again,  the  first  pulse  in  the  train  excites  carriers  to  the  excited 
state  absorption  level  but  the  carriers  do  not  have  sufficient  time  to  return  to  the  ground 
state  before  the  next  pulse  arrives.  This  leads  to  nearly  constant  transmission  energy  for 
the  trailing  pulses. 

We  have  numerically  investigated  the  effects  of  multiple  pulses  on  excited-state 
optical  limiters.  For  closely  spaced  pulses,  pulse  interactions  leads  to  enhanced  peak 
energy  and  improved  limiting  behavior.  However,  even  for  widely  spaced  pulses,  the 
trailing  pulses  reach  a  near  constant,  but  lower,  transmission  energy.  Our  initial  results 
show  that  transmission  of  pulse  trains  behave  differently  than  transmission  of  a  single 
pulse.  Since  many  lasers  produce  pulse  trains  or  are  externally  modulated  to  produce 
pulse  trains,  the  effects  of  multiple  pulses  on  optical  limiting  are  of  increased  importance. 

7.0  Conclusions 

Increasing  usage  of  high  power  lasers  in  laboratoiy  research,  measurement,  medical 
treatment  and  weapons  make  human  eyes  and  other  optical  imagining/detecting  systems 
more  vulnerable  to  laser  related  damage.  Laser  induced  detector  damage  such  as  electron 
breakdown  in  semiconductor  detectors  or  condensed  media,  thermal  acoustic  explosion 
and  thermal  chemical  reaction  in  the  retina  can  happen  when  the  detector  is  illuminated 
by  a  very  high  irradiance  laser  beam  directly  or  indirectly  (through  some  kind  focusing 
optical  structure).  An  effective  optical  limiting  device  can  prevent  or  reduce  damage  by 
clamping  the  excessive  energy  and  power  from  these  high  intensity  light  sources. 

Nonlinear  properties  including  nonlinear  refractive,  absorption  and  scattering  make 
passive  optical  limiting  possible.  Over  the  last  several  decades,  many  materials  have  been 
discovered  that  exhibit  optical  limiting  properties.  These  materials  include 
semiconductors  (TP A),  carbon-black  suspension  (thermal  scattering),  fullerenes  (RSA), 
and  porphyrins  (RSA).  Currently  porphyrins  and  their  derivatives  are  considered  by  many 
to  be  very  promising  materials. 

The  nonlinear  mechanisms  used  for  optical  limiting  have  different  properties  that 
require  various  models.  For  example,  nonlinear  scattering  in  CBS  is  fluence  dependent 
with  broadband  features,  Kerr  nonlinearity  and  two-photon  absorption  are  intensity 
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dependent  and  occur  instantaneously,  while  the  nonlinear  absorption  in  RS  A  material  can 
be  fluence  or  intensity  dependent  and  depends  on  pulse  width  and  beam  shape. 

One  of  the  key  features  to  the  successful  design  of  optical  limiting  is  numerical 
simulation.  The  increasing  complexity  of  the  materials,  device  geometry,  and  optical 
beam  shaping  and  signal  processing  necessitates  the  use  of  numerical  simulation  for 
system  optimization.  In  an  optical  limiting  system  using  RSA  material,  the  propagation 
equations  are  coupled  to  the  rate  equation.  The  mathematical  methods  used  for  optical 
limiter  calculations  include  the  Huygen-Frensel  Principle,  the  Crank-Nicholson  finite 
difference  method,  the  spectral  or  Fourier  transform  method  and  the  split-step  method. 

The  effects  of  spatiotemporal  behavior  in  optical  limiting  systems  have  only  recently 
been  addressed.  However,  these  phenomena  could  potentially  play  a  significant  role  in 
optical  limiting.  For  example,  in  the  eye,  pulse  compression  due  to  nonlinear  Kerr  effect 
and  negative  GVD  may  reduce  the  damage  threshold;  whereas,  the  Kerr  effect  and 
positive  GVD  may  increase  the  damage  threshold  in  the  subpicosecond  region.  Until 
recently,  calculations  have  not  included  the  effects  of  both  nonlinear  Kerr  effect  and 
dispersion.  The  nonlinear  Kerr  effect  and  diffraction  may  have  significant  effects  on  the 
optical  limiting  system.  Self-focusing  can  change  the  focal  point  in  the  system. 
Additionally  self-focusing  can  shift  the  focal  point  to  the  exit  plane  of  the  nonlinear 
absorbing  material  and  increase  its  optical  limiting  property  by  utilizing  all  of  the 
molecules  in  the  optical  path. 

Though  some  optical  limiters  have  already  been  used  in  real  applications,  they  are  still 
insufficient  to  protect  very  sensitive  devices,  like  human  eyes,  from  very  high-energy 
beams.  New  materials  with  high  nonlinearity  and  fast  response  will  remain  a  key  research 
area.  Additionally  investigations  on  the  interplay  and  synergism  between  different  types 
of  nonlinear  optical  phenomena  may  significantly  improve  the  performance  of  optical 
limiting  systems.  New  applications  of  these  optical  limiting  materials  will  also  be  an 
interesting  area. 
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Fig.l.  Schematic  of  approximate  time  scales  for  various  optical  limiting  phenomena. 
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Fig.  2.  Schematic  energy  level  for  a  chromophore.  The  electronic  states  are  represented 
by  solid  horizontal  lines  and  the  vibronic  states  are  represented  by  dotted  horizontal  lines. 
Sq  represents  a  singlet  state  and  Tq  represents  a  triplet  state.  The  photon  absorption 
excitations  are  represented  by  solid  vertical  lines,  and  the  decay  processes  are  represented 
by  wavy  lines.  The  absorption  coefficients  apq  and  the  decay  constants  kpq  are  described 
in  the  text.  The  physical  values  used  for  our  calculations  are:  ctoi=2.4x10'18  cm2, 
(712=3.0x1  O'17  cm2,  <734=4.8x10'17  cm2,  kio=0.144  ns*1,  k2i=1.0  ps*1,  ki3=77.8  pis'1, 
k3o=50.0  ms'1, 1(43=1.0  ps'1,  and  kn,  k22,  k33  are  due  to  vibrational  decays  which  are 
assumed  to  be  instantaneous. 
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Fig.  3 .  Diagram  of  the  features  of  the  numerical  beam  propagation  code. 
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Fig.  4.  Pulse  shapes  of  the  laser  beam  at  four  propagation  distances.  The  incident  energy 
is  10  pJ  and  the  carrier  density  Nt=6.022x1017  cm'3. 
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Fig.  5.  Carrier  densities  as  a  function  of  space  and  time,  (a)  Ground  state:  Maximum 
value  =1.0  and  minimum  value  =0.772;  (b)  First  excited  state  of  single  state:  Maximum 
value=0.226;  (c)  Second  excited  state  of  single  state:  Maximum  value=0. 00624;  (d)  First 
state  of  the  triplet  state:  Maximum  value=0.00466. 
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Peak=0.0 14,  Total=0.01 8  Peak=0.002,  Total=0.007 


Fig.  6.  Pulse  shapes  of  the  laser  beam  at  four  propagation  distances.  The  incident  energy 
is  500  pJ  and  the  carrier  density  Nt=6.022x1017  cm'3. 


Incident  Energy=500  \jJ,  Z=0.28  cm 
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Fig.  7.  The  carrier  densities  at  z=0.28  cm  and  incident  energy=500  |jJ. 
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Fig.  8.  Intensity  and  carrier  density  plots  for  input  energy  =  50  joJ:  (I.a)  the  dotted  line 
shows  the  normalized  input  pulse,  the  solid  line  shows  the  output  pulse  whose  peak  value 
=0.40;  (Lb)  the  dashed  line  shows  the  carrier  density  of  state  So,  whose  minimum  value  is 
0.13,  the  solid  line  shows  the  carrier  density  of  state  Si  whose  maximum  value  is  0.86, 
the  dot-dashed  line  shows  the  carrier  density  of  S2  whose  maximum  value  is  0. 19;  (jQ) 
shows  the  3-D  plots  of  the  carrier  densities  at  (a)  the  ground  state,  (b)  the  first  excited 
singlet  state,  (c)  the  second  excited  singlet  state,  (d)  the  first  excited  triplet  state. 
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Fig.  9.  Intensity  and  carrier  density  plots  for  input  energy  =  1000  pJ:  (La)  the  dotted  line 
shows  the  normalized  input  pulse,  the  solid  line  shows  the  output  pulse  whose  peak  value 
=0.72;  (Lb)  the  dashed  line  shows  the  carrier  density  of  state  So  whose  minimum  value  is 
0.0,  the  solid  line  shows  the  carrier  density  of  state  Si  whose  maximum  value  is  0.99,  the 
dot-dashed  line  shows  the  carrier  density  of  S2  whose  maximum  value  is  0.99;  (II) 
shows  the  3-D  plots  of  the  carrier  densities  at  (a)  the  ground  state,  (b)  the  first  excited 
singlet  state,  (c)  the  second  excited  singlet  state,  (d)  the  first  excited  triplet  state. 


Fig.  10.  Transmission  as  a  function  of  the  coefficient  of  the  nonlinear  index  of  refraction 
(n2)  in  units  of  10'17  cmVfiW.  The  incident  energy  is  1  |jJ  and  the  carrier  density  is 
Nt=6.022x1016  cm'3.  The  increased  absorption  is  due  to  the  self-focusing  of  the  light 
arising  from  the  nonlinear  Kerr-effect. 
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Figure  11.  (a)  Kerr  self-focusing  (b)  tandem  optical  limiter  using  self-focusing  before 
nonlinear  absorber. 


Figure  12.  Schematic  diagram  of  transmission  of  tandem  optical  limiter  using  self- 
focusing. 
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Figure  13.  Pulse  reshaping  in  a  tandem  optical  limiter  using  Kerr  nonlinearity  and  GVD 
in  element  (1)  and  nonlinear  absorption  in  element  (2);  (a)  pulse  compression,  (b)  pulse 
broadening. 
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Fig  14.  (a)  Gaussian  input  pulse  with  its  peak  value  normalized  to  1.  (b)  The  pulse  shape 
after  propagating  one  diffraction  length  in  a  nonlinear  dispersive  material  (NL=2.0, 
GVD--0.05);  its  peak  value  is  18  even  though  in  the  graph  z  axis  value  is  normalized  to 
1. 
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Fig.  15.  Carrier  density  population  profiles  at  a  propagation  distance  of  one  diffraction 
length  in  the  nonlinear  material  and  0.7  mm  in  the  RSA  material.  The  input  energy  is  50 
uJ.  (a)  The  ground  level,  min.  value=  2.3e-5.  (b)  The  first  excited  singlet  (SI)  state,  max. 
value=  0.984  (c)  The  second  excited  singlet(S2)  state,  max  value=0.912.  (d)  The  first 
excited  triplet(Tl)  state,  max  value=  0.019. 
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Fig.  16.  Input  energy  vs.  transmission  for  of  one-diffraction  length  (NL  +GVD)  and  0.72 
mm  propagation  through  RSA.  The  circle  data  points  are  for  NL=0.0,  GVD=0.0  case. 
The  square  data  points  are  for  NL=2.0,  GVD=-0.05  case. 
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Fig.  17.  Schematic  diagrams  of  possible  optical  beam  processing  for  numerical 
calculation  of  optical  limiting:  (a)  multiple  input  pulses;  (b)  randomization  of  temporal 
pulses  using  a  modulator. 
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Fig.  18.  Transmission  versus  input  energy  for  xr  =2  and  10. 
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Fig.  19.  Normalized  intensity  curves  for  input  energy  =  9|_i  J  and  tr=2;  (a)  and  (b)  before 
excited-state  absorber;  (c)  and  (d)  after  excited-state  absorber. 
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Fig.  20.  Normalized  intensity  curves  for  input  energy  =  9\i  J  and  xr=10;  (a)  and  (b)  before 
excited-state  absorber;  (c)  and  (d)  after  excited-state  absorber. 
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